Interior gap superfluidity in a two-component Fermi gas of atoms 
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A new superfluid phase in Fermi matter, termed as "interior gap" (IG) or "breached pair", has 
been recently predicted by Liu and Wilczek [Phys.Rev.Lett. 90, 047002 (2003)]. This resuhs 
from pairing between fermions of two species having essentially different Fermi surfaces. Using 
a nonperturbative variational approach, we analyze the features, such as energy gap, momentum 
distributions, and elementary excitations associated with the predicted phase. We discuss possible 
realization of this phase in two-component Fermi gases in an optical trap. 
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Since the first realization of Fermi degeneracy in K 
gas by Jin's group 1| , several other groups "5, 3, P 



have also achieved quantum degeneracy in atomic Fermi 
gases. A recent experiment by Regal et at 7], has re- 
ported condensation of fermionic atom-pairs in trapped 
""^K gas. Recently, several groups have demonstrated 
formation of cold molecules [3, |9( and molecular Bose- 
Einstein condensates from trapped Fermi gases of 
atoms near the Feshbach resonance. O'Hara et al. 
demonstrated strongly interacting degenerate Fermi gas 
of ^Li atoms. Studies on trapped degenerate Fermi gases 
are important in the context of diverse fields, such as 
nuclear physics, astrophysics, strong-coupling supercon- 
ductivity and superfluidity, and so on. 

A new form of Fermi matter, associated with what 
is termed as "interior gap " (IG) superfluidity llSj , has 
recently come into fore. Liu and Wilczek [fj] have 
shown that, in two species of fermions differing in Fermi 
momenta, an attractive interaction can lead to pairing 
within the interior of the larger Fermi surface, while the 
exterior remains gapless. This pairing phenomenon gives 
rise to IG superfluidity having superfluid and normal 
Fermi liquid components simultaneously. Many years 
ago, an analogous state of two-species fermions was pre- 
dicted by Sarma |l4j . which, however, corresponded to a 
metastable phase. Recently, it was shown in the context 
of color superconductivity that, such a phase is the only 
sta ble p hase when the relative density of the two species 
[isl is kept fixed. IG superfluidity arises naturally in 
finite density quark matter [11, lH, ■ 

In this paper, we analyze the energy gap, momen- 
tum distributions and quasiparticle excitations in an 
IG superfluid system. For possible realizations of this 
phase, optically trapped two-component atomic Fermi 
gases may be useful. It has already been suggested 
that a mixture of ^Li and "^^K Fermi gases should be used 



to create this new phase. Let us consider a mixture of 
two hyperfine spin components, for instance, |1) = |_F 
1/2, mp = 1/2) and |2) = |F = 1/2, = -1/2) of *^Li 
atoms, in an optical trap. The relative number density 
of the two components can be controlled by a rf field 
|^ii,2IL 21|. The lifetime of each spin component is long 
enough to carry out an experiment at a fixed relative 
density which is a necessary condition for IG superfluid- 
ity [131 ■ It may be noted that controlled mixtures, other 
than 50:50, of fermion gas have already been created [2l| . 

To discuss IG superfluidity, let us consider a system of 
interacting fermions of two species in a harmonic trap, 
described by the Hamiltonian H = Hq + Hj where 
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where Vhoij) is the harmonic trapping potential and 



Hr = 



drdr'#t(r)'i'^(r')F(|r 
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r'|)*i(r)vl/2(r'). 

(2) 

Here, 5'i(r) represents fermionic held operator, /ii(r) 
denotes the local chemical potential and T^(|r — r'|) 
is the two-body interaction potential. We now in- 
troduce the annihilation operator Ci(k) which, for a 
uniform gas, is related to field operator by 'l'i(r) — 
{l/\/V)^^exp{—ik.r)ci{k), where V is the volume of 
the system. For a uniform gas with number density n, 
the mean field energy in the weak interaction regime 
is Umf = gn, where g = Airh'^as/ (2rh), with rh = 
mim2/{mi + being the reduced mass, and as the s- 
wave scattering length. In the strong interaction regime, 
it is unitarity-limited and takes the form Umf — (3Ep 
[ill I23I li^. where is a constant. In this regime, one 
can define an effective scattering length from the rela- 
tions Umf — gn = PEp yielding UoS oc 1/kp. Under 
local density approximation (LDA) in a harmonic trap 
of potential Vho, the equilibrium conditions are given by 
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— [67rV(r)]2/3+g|e,,>,(r) + y,o(r) =/i,(r), (3) 

where is the Levi-Civita tensor. LDA is particularly 
good when the local Fermi energy is larger than the av- 
erage level spacing of the trap and the coherence length 
of the fermion pair is shorter than the average trap size 
[2^ . The Fermi momentum of the ith component is then 
related to the respective peak density ni(ro) at trap cen- 
ter by kp. = [6n'^ni{ro)]^^^ ■ In what follows, we thus 
consider a uniform gas and derive the gap equation. 

To consider the pairing state, we take the variational 
ansatz as given by |l8l l26l | 

\n)^exp (c,(k)t/(k)c,(-k)t)e,,dk-H.c. |0), 

where |0) represents the vacuum, annihilated by Ci's. 
To include the effects of temperature and density, we 
write down the state at finite temperature and density 
|r2(/3,^)) taking a thermal Bogoliubov transformation 
over the state |r2) using thermofield dynamics as de- 
scribed in Refs. HlUll. We then have, 

=exp(^y c^(k)t0,(k,/3)4(-k)tdk-H.c.^ |r!>. 

(5) 

The superscript "prime" refers to operators in basis, 
i.e., c^\ft) = 0. The tilde operators are the ones in the ex- 
tended Hilbert space. In Eq. ((SJ, the function 6'i(k, f3, fj) 
will be related to the distribution functions, for the two 
species. Using the Bogoliubov transformation, it is easy 
to calculate the thermodynamic potential given by 



(6) 



where Ti. is the Hamiltonian density for con- 
stant potential. In the above [23|, s = 

- J2tl d'k(^ sin^ 9, In sin^ 9^ + cos^ 9i In cos^ 9i ^ 
is the entropy density. 

Minimizing the thermodynamic potential f2 with re- 
spect to condensate function /(k), we get tan2/(k) = 
-{2glD)/{ei - Ml + gp2 + £2 - M2 + 9Pi) = A/(e - v), 
where we have defined the gap A = —gljj. The order 
parameter here is (V'i(x)-02(x)) = —Id where, Id — 
/ sin2/(k) (1 - sin^ 9i - sin^ 92) dk. Further, e = 
(ei -|- £2)72, V = {vi + i^2)/2, Vi is the interacting chem- 



ical potential given as Vi 



Pi 



9Pj\ 



Thus, it 



may be noted that the pair condensate functions de- 
pend upon the average energy and the average chem- 
ical potential of the fermions that condense. Finally, 
the minimization of the thermodynamic potential with 
respect to the thermal functions 9i (k) yields sin'^ 9i = 
l/(exp(/3a;i) -I- 1). The quasi particle energies w^'s are 



given hy LUi — oj + S^ L02 = — Here, oj = \/ -I- f ^ , 
^ = e— V and (5^ = (ei — e2)/2 — [vi — 1^2)/2. Note that in 
the degenerate mass and chemical potential case we shall 
have the same quasi particle energies but if the masses 
and/or the chemical potentials are different there is pos- 
sibility of having gapless modes wi = {uj2 = 0) when 
Lj = —6^ ( uj ~ d^). So, although we shall have nonzero 
order parameter A, there will be fermionic zero modes or 
gapless modes [iTlllg- Indeed, we shall show below that 
such a situation is possible depending upon the mismatch 
in the fermion number densities and the magnitude of A. 

Substituting the solutions for the ansatz functions in 
definition for Id, we can have the gap equation which, 
however, is divergent. The origin of this divergence lies 
in the point-like four fermion interaction which needs to 
be regularized. We tackle this problem by defining the 
regularized coupling as was done in Refs. j23,|22| to access 
the strong coupling regime 29, 30] by subtracting out 
the zero temperature and zero density contribution. The 
regularized gap equation is then given by 
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{2n) 



dk 



2VA2 -f ^2 2e 



(7) 



where gu = 2iTh^aes/'fn. Unless otherwise stated explic- 
itly, we take /i = w^? = 1 in our calculations, where vp 
is the Fermi velocity. Substituting the expressions for 
the optimized condensate functions and the distribution 
functions, finite part of the thermodynamic potential Eq. 
in the zero temperature limit, reduces to 



1 



(27r)3 



dk 



A2 

9R 



9B.P1P2 



dk [{lo + 8e)9{-LOi) + {uj- 6(_)9{-uj2)] (8) 



In the case of degenerate masses and chemical poten- 
tials, the theta functions will not contribute. The solu- 
tion of Eq. (0 is shown in Fig. In the dilute gas limit 
(fc_F|acff| « 1), we have the BCS solution |2^Eg| An c± 
8i?_Fe"2exp[-7r/(2fco|aoff|)], where Ep = kl/Arh. This 
can be seen in Fig. 2(a) where A/Ep rises exponentially 
for small fco|ae//|. At a very high value of fcn|Qf,ff|, the 
gap becomes proportional to the Fermi energy [23 . l29l| . 
For the case when is nonzero, one of the two theta 
functions in Eq. {T)) will contribute. To be specific let us 
consider the case when 6^ is negative. In that case the 
contribution from 9{—lui) will be nonzero in a strip of mo- 



menta defined by fc^in < k < k^. 
{miVi^'m2V2)±\/{miVi - m2V2Y 
The gap equation now becomes 



where fc^ax.mm 



1 = 



\9R\ 
(27r)3 
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v/a^TF 
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dk, 
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FIG. 1: Right hand side (RHS) of gap Eq. is plotted vs 
A|oeff| for different values of difference in chemical potential 
5^. The gap is given by the condition RHS=1. The solid, dot- 
ted dashed-dotted and dashed lines correspond to (5^|acff |=0, 
0.25, 0.35 and 0.45, respectively. For (5^|acff| ~ 0.45, two 
solutions for gap exists. The smaller value corresponds to 
"IG state" and the larger one corresponds to usual Bardeen- 
Cooper-Schrieffer state. 



where Skp — kp 



For fixed average chemical poten- 
tial and for nonzero 5^ there are, in general, two solutions 
for the gap as shown in Fig. 1. Clearly, when ^ 0, 
there is only one solution which for the parameters chosen 
turns out to be A = 0.71\armef f\~^ ■ As is increased, 
the second integral in Eq. © starts contributing and 
the curve can cross in general the RHS=1 line at two 
places. The larger A is the conventional BCS solution. 
The smaller value of A will have breached gap character 
(isf and, in this case, the number densities of the two 
species will be different. In fact, the number densities of 
the two species are given as, for 5^ < 0, 



1 



Pi = 



(27r)3 



P2 = 



(27r)3 



1 - 



{1 - e{-Lui)) dk. 



dk, 
(10) 

(11) 



Clearly, without any mismatch in the chemical poten- 
tial of the two species, loi — lu > and hence the two 
densities are equal as the theta functions will not con- 
tributes and the distribution function becomes the usual 
BCS distribution function. In presence of mismatch in 
densities, the 9 function contributes in the range of mo- 
menta between kmin and kmax- In that region, the parti- 
cle distribution function becomes unity for species one 
and zero for the other species. This is what is seen 
in Figs. 2(c) and 2(d). It is easy to see that in the 
breached gap phase, the difference in number densities is 



FIG. 2; (a) The gap (in unit of Ef) as a function of fco|acff|. 
(b) The gap as a function of ratio [p\/ P2) of the two densities. 
Average occupation numbers ni(fc) and n2{k) as functions of 
k are shown in (c) and (d), respectively, for the parameters 
fcolfflcff I = 1.5 and pi/p2 = 1.4. 



For systems where chemical potential is kept fixed, 
the breached pair phase will always have higher free 
energy as compared to the BCS phase. However, for 
systems where the relative densities of the condensing 
fermions can be fixed, the breached pairing state will be 
the only solution provided the ratio between the densi- 
ties is different from unity. In Fig 2(b), almost lin- 
ear behavior of the gap with the density ratio may be 
noted. Figures 2(c) and (d) show that in a strip of mo- 
menta between kmin and kmax, the momentum distri- 
bution resembles that of normal fermionic distribution. 
Below kmin and beyond kmax, the distribution is of the 
BCS type. The quasi-particle excitation energy u>i be- 
comes zero at the points /cmin and fc,„ax showing gap- 
less behavior. For the parameters chosen in Figs. 2(c) 
and (d), A:i„in|aeff| = 1.42 and /cmaxlacffi = 1-58, i.e., 
kmin/ko = 0.95 and k^i^^/ko — 1.05, fco being the av- 
erage Fermi momentum of the two species. 

We thus show that the ansatz state, as defined in Eq. 
(jsj, shows breached pair superfluidity with the excita- 
tion energies becoming zero at k^i^ and fcmax and that 
when relative number density is fixed this state is the 
state having lower free energy. In this region, only one of 
the species is completely occupied and the other is com- 
pletely empty. Thus, it is a more general state as com- 
pared to Ref.^3 where the excitation energy was zero at 
/cniin which was equal to pA of Ref. beyond which 

one species was completely occupied and the other was 
completely empty. 

We next discuss the possibility of probing the gap by 
stimulated Bragg scattering (SBS) 31, 32^ in analogy 
with Raman scattering in anisotropic superconductors 
[33! . There are many theoretical proposals 0, for 
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FIG. 3: A possible scheme for measuring dynamic structure 
function by SBS. The laser photon with momentum ki is 
scattered into a photon with momentum ks transferring a 
momentum q = ki ~ ka to the atoms along the x axis, while 
the laser photon with momentum ka is scattered into a photon 
with momentum k2 transferring a momentum — q = ka — k2 to 
the atoms. The frequencies of the three laser beams are chosen 
such that the energy transfer 5 = (oji — ws) = (wa — tJ2) > 0. A 
pair of atoms which may form a Cooper pair having mutually 
opposite momenta are scattered equally. The one-dimensional 
(along the x direction) momentum and density distribution of 
the gas can be determined from the analysis of time of flight 
images. The spectrum (number of scattered atoms versus S 
for different q values) may reveal the existence of the gap and 
the detailed comparison of the spectra for different relative 
densities of the two components may provide a proof of IG 
superfluidity. 



probing atomic Cooper pairs by resonant |2J| and ofF- 
resonant js^ laser light. It is well established that, in 
the case of a Fermi superfiuid, there will be no spectrum 
of density fluctuation unless energy transferred to the 
superfiuid due to scattering exceeds 2 A [s^- The quan- 
tum Monte Carlo simulation suggests that, in the 
large limit, A ~ 0.81(3/5)e_F — 0.49eF- In the recent 
experiments [3, Hol llH on two-component ^Li atoms in 
the unitarity- limited regime (|as|fci? » 1), the typical 
value of the Fermi velocity vp = hkp /m ~ 15 cm/s. For 
two-component ^Li atoms, if the Bragg pulses are tuned 
near an excited level (say, 2P3/2, wavelength A — 670.776 
nm), the momentum transfer q = ki — ka ~ ik^'k, where 
fci = 2'k/X. This momentum transfer raises the velocity 
of the scattered atom by 2(fifci/m) ~ 20 cm/sec which 
exceeds vp- Therefore, the scattered atoms can be dis- 
tinguished in time of flight images as discussed in the 
caption of Fig. 13 

In conclusion, we have analyzed momentum distri- 
butions and energy gap for an IG superfiuid. Two- 
component Fermi gases of atoms in optical traps seem 
to be promising systems for experimental exploration for 
such superfiuid state. Fermi atoms in a two-band optical 
lattice may also be considered for studying breached pair 
superfluidity. 
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